Enhanced Dimer Relaxation in an Atomic/Molecular BEC 
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Using Feshbach resonances it is possible to study 
cold atoms with a variable scattering length a that can 
be made large enough that the atoms have strong reso- 
nant interactions. Surprising phenomena have been ob- 
served in experiments on cold atoms with large scatter- 
ing lengths. Experiments with a 23 Na BEC 0, ultracold 
85 Rb atoms 0, a 87 Rb BEC Q, and ultracold 133 Cs 
atoms [B| have revealed large losses of atoms as the Fes- 
hbach resonance is approached. In experiments with a 
85 Rb BEC in which a was increased briefly to large posi- 
tive values, a burst of energetic atoms was observed along 
with missing atoms 0. More recently, atom-molecule 
coherence was demonstrated [7], suggesting the coexis- 
tence of a condensate of shallow diatomic molecules with 
an atom condensate. Few-body processes play an im- 
portant role in these phenomena and it is important to 
understand their behavior in the limit of large a. 

A Feshbach resonance allows \a\ to be made large 
compared to the typical low-energy length scale I = 
(mCe/h 2 ) 1 ^ 4 set by the atom mass m and the van der 
Waals coefficient C§. If \a\ ^> £, the few-body problem 
exhibits universal properties that are insensitive to the 
details of the interactions responsible for the large scat- 
tering length. A simple example in the 2-body sector 
for a > is the existence of a shallow molecule (the 
dimer) with binding energy \E 2 \ = h 2 /{ma 2 ). A par- 
ticularly remarkable example in the 3-body sector is the 
existence of shallow 3-bodybound states (Efimov states) 
with universal properties 8J. As \a\ increases, shallower 
Efimov states appear below the scattering threshold at 
values of a that differ by multiples of e 71 ^ rs 22.7, where 
s sa 1.00624. The Efimov states can cause dramatic 
dependence of scattering observables on a and on the en- 
ergy. The 3-body scattering observables are determined 



by a and by a 3-body parameter that also determines 
the Efimov spectrum. Simple universal expressions have 
been derived for the atom- dimer scattering length 0, |j| 
and for the rate constant for 3-body recombination into 
the shallow dimer [ToL ITU Il2^ . They are valid for all 
atoms with \a\ ^> I that have no 2-body bound states if 
a < and only the shallow dimer if a > 0. 

Unfortunately, heavy alkali atoms, such as Rb and Na, 
form many deeply-bound diatomic molecules. A simple 
consequence is that Efimov states are resonances rather 
than sharp states, because they can decay into a deep 
molecule and an energetic atom. The widths of Efimov 
resonances have been studied in Refs. using explicit 
model potentials. The deep molecules also affect other 3- 
body observables. In particular, the universal expressions 
for the atom-dimer scattering length and for the 3-body 
recombination rate constant must be modified. 

The existence of deep molecules opens up additional 
channels for the loss of cold atoms from an ultracold gas 
or a BEC. One such channel is 3-body recombination 
into dee p m olecules, which has been studied previously 
in Refs. [ill [lj] • Another loss mechanism for a cold gas 
containing both atoms and dimers is dimer relaxation, 
in which the collision of an atom and a shallow dimer 
produces a deep molecule and a recoiling atom. If the 
momenta of the incoming atom and dimer are sufficiently 
small, the relaxation rate reduces to a constant. The 
relaxation event rate constant [3 is defined so that the 
number of relaxation events per time and per volume in 
a cold gas (or BEC) of atoms and dimers with number 
densities ua and no is PnAno- Because of the large 
energy released, this process results in the loss of a low- 
energy atom and a shallow dimer: 
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dt dt 

The loss of atoms from a BEC through dimer relaxation 
has been considered previously in Refs. 0] . The authors 
assumed that the relaxation rate constant [3 is indepen- 
dent of a. They found qualitative agreement with the 
slow-sweep experiments in 2] for j3 ~ 10~ 10 cm 3 /s. 
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In this paper, we derive a universal expression for the 
dimer relaxation rate constant (3 in the limit of large 
a using hyperspherical methods combined with previous 
results from an effective field theory. We show that j3 
scales like ha/m, but with a coefficient that is a periodic 
function of ln(a) with period tt/sq ln(22.7). The coef- 
ficient of Ha/m also depends on two 3-body parameters 
that determine the binding energies and widths of the 
Efimov resonances. It can have resonant enhancement at 
those values of a for which there is an Efimov resonance 
near the atom-dimer threshold. 

In contrast to other treatments of Feshbach resonances, 
we do not attempt to describe the details of the physics in 
the short-distance region. Such details include the closed 
channel of the Feshbach resonance and the decomposition 
of the dimer relaxation rate into contributions from indi- 
vidual deep 2-body bound states. A detailed description 
of the short-distance physics is not necessary for calcu- 
lating the inclusive dimer relaxation rate. Since the total 
probability is conserved, one can describe the inclusive 
rate by taking advantage of unitarity. 

We first consider the case without deep dimers. Any 
probability that flows to the short-distance region where 
the atoms are coupled to the closed channel must be 
reflected back to the long-distance region in the form 
of low-momentum states. A detailed knowledge of the 
short-distance dynamics of the Feshbach resonance which 
creates the large scattering length is not required. One 
can eliminate the closed channel in favor of the scatter- 
ing length a and the 3-body parameter A* discussed be- 
low. This simplification is possible because the resolution 
scale of the low-momentum atoms is given by the scat- 
tering length a which is much larger than the range of 
the potential for the closed channel. As a consequence, 
low-energy atoms are only able to feel the large scatter- 
ing length a (and the three-body parameter A* if three- 
and higher-body observables are concerned) but not to 
probe the mechanism responsible for it. 

If deep dimers are present, some of the probability that 
flows to the short-distance region will emerge in the form 
of high-momentum states consisting of an atom and a 
deep dimer. However, if one is interested in the inclusive 
relaxation rate only, it is irrelevant how exactly the prob- 
ability is lost from the low-momentum states. All one 
needs to know is how much of the probability is lost. In 
the case of atoms with a large scattering length, this can 
be parametrized by one number 77*. The reason for this 
is again that the resolution scale of the low-momentum 
atoms is given by the scattering length a which is much 
larger than the size of the deep dimer. Therefore the 
low-momentum atoms cannot resolve the details of the 
physics at short distances comparable to the size of the 
deep dimer. 

A similar approach was previously used to describe 
the effect of deep 2-body bound states on Efimov states 
fl(| or the three-body recombination rate [l4|- It has 
also been applied to many other problems in physics. In 
the mean- field equations for Bose-Einstein condensates, 



atom loss processes are often taken into account through 
imaginary parts of coefficients in the mean-field Hamil- 
tonian 17]. If the loss processes involve momenta large 
compared to those of atoms in the condensate, this pre- 
scription can be rigorously justified. A similar method 
is used in nuclear and particle physics to describe ef- 
fects of the annihilation of particles and their antiparti- 
cles 01 • The dynamics of the short-distance annihilation 
process may be very complicated, but its effects on low- 
momentum particles can be described by local operators 
in an effective Hamiltonian and all one needs to know is 
the imaginary parts of the coefficients of those terms. 

In the remainder of the paper, we will make these ar- 
guments more quantitative and use them to derive a uni- 
versal expression for the dimer relaxation rate constant (5 
in the limit of large a. We begin by considering the equa- 
tion for the hyperradial wavefunction f(R) for the lowest 
adiabatic potential in the hyperspherical representation 
of the 3-body problem [Tjj. In the region i -C R -C a, 
the other adiabatic potentials are repulsive and the corre- 
sponding wavefunctions decrease exponentially at R <C a. 
The coupling to those channels can therefore be neglected 
and the equation for f(R) reduces to 
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f(R) = Ef(R) 



(2) 



This looks like the Schrodingcr equation for a particle in a 
one-dimensional, scale- invariant 1/R 2 potential. For low 
energies \E\ ~ ft 2 /ma 2 , the energy eigenvalue in Eq. 
can be neglected and the most general solution is 



f(R) 



(3) 



where H has dimensions of wave number but is otherwise 
arbitrary. Equation J3J) is simply the sum of outgoing and 
incoming hyperradial waves. At shorter distances R ~ I 
and longer distances R ~ a, the wave function becomes 
very complicated, but it is constrained by unitarity. 

Following Efimov we first exploit unitarity in the 
long-distance region R ~ a. If we identify the asymp- 
totic states at distances R <C a and R^> a, the evolution 
of the wave function through the region R ~ a can be 
described by a unitary matrix s (see the illustration in 
Fig. Pi . We denote the asymptotic states with probabil- 
ity flowing into and out of the region R ~ a by kets \i in) 
and \i out). The probability can flow into this region 
either from £ <C R <C a or from R ^> a. In the region 
< C fi € a, the asymptotic states |1 in) and |1 out) are 
the outgoing and incoming hyperradial waves represented 
by the first and second terms in Eq. (jJJ. For simplicity, 
we restrict ourselves to energies below the dimer-breakup 
threshold where 3-atom states are kinematically forbid- 
den. In the region R ^> a, the asymptotic states |2 in) 
and 1 2 out) are incoming and outgoing atom-dimer scat- 
tering states. The amplitudes for incoming asymptotic 
states \j in) to evolve into outgoing asymptotic states 
\i out) form a unitary 2x2 matrix . Time-reversal in- 
variance implies that s is a symmetric matrix. Since s is 
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FIG. 1: The hyperspherical potential in the region I <C R <C 
a. The evolution of wavefunctions in the regions R ~ a and 
R ~ Z is described by unitary matrices s and i. 

dimensionless, it can depend on the interaction potential 
and the energy only through the variable ma 2 E/h 2 . 

Following Efimov, we next exploit unitarity in the 
short-distance region R ~ £. If there are no deep 
molecules, the probability in the incoming hyperradial 
wave must be totally reflected from the region R ~ £ 
0. Thus the amplitudes of the incoming and outgoing 
hyperradial waves in J2J must be equal in magnitude: 
A = —Be 2ie * . The angle is determined by the loga- 
rithmic derivative f'(R)/f(R) in the region ii«a. This 
implies that 6* must have the form 8* = so ln(CaA*), 
where C is a dimensionless numerical constant and A* is 
a parameter that is insensitive to variations in a. The 
evolution of the wave function in the region R ~ £ can 
be described by a 1 x 1 unitary matrix whose only entry 
is tn = e 2lS * . By combining the constraints from uni- 
tarity in the regions R ~ £ and R ~ a, Efimov obtained 
an expression for the S-matrix element for atom-dimer 
scattering that he referred to as a radial law |8|]: 

SAD,AD = S 2 2 + «2l(l - *llSl2 . (4) 

Because s is dimensionless, it can only depend on the di- 
mensionless variable ma 2 E/h 2 . The atom-dimer scatter- 
ing phase shift can therefore be obtained by calculating 
a few universal functions of ma 2 E/Ti 2 

We now consider the case with deeply-bound diatomic 
molecules. In this case, some of the probability in an in- 
coming hyperradial wave that flows into the region R ~ £ 
will emerge as high-energy scattering states of an atom 
and a deep molecule. The amplitude of the outgoing hy- 
perradial wave in © must therefore be smaller than that 
of the incoming hyperradial wave: A = —Be 2 ' te '~ 2v *. 
The corresponding element of the unitary matrix t is 

tn = exp(2j6>» — 2?7*) . (5) 

Inserting this expression into (£Q), we obtain the general- 
ization of Efimov's radial law to the case in which there 
are deep molecules. If the universal expressions for the 



S-matrix elements are known in the case 77* = 0, the 
effects of deep molecules can be deduced by the simple 
substitution 6>* — > 9* + irj* . 

We proceed to derive the radial law for S-matrix ele- 
ments involving the deep molecules. In the short-distance 
region R ~ £, the relevant asymptotic states are incom- 
ing and outgoing hyperradial waves |1 in) and |1 out) 
and high-energy atom-molecule scattering states which 
we denote by \X in) and \X out). Note that the out- 
going hyperradial wave is an incoming asymptotic state 
|1 in) as far as the region R ~ a is concerned and an 
outgoing asymptotic state |1 out) as far as the region 
R ~ £ is concerned. The S-matrix element for a low- 
energy atom-dimer scattering state to evolve into a high- 
energy atom-molecule scattering state X is 

Sx,AD = t X \Sl2 + t X lSu(l - <llSll) _1 illSi2 . (6) 

The first term on the right side of Eq. © describes trans- 
mission through both the regions £ <C R <C a and R ~ £ 
to the asymptotic state X . If we expand the second term 
in powers of tuSn, we can identify the n th term as the 
contribution from transmission through the region R ~ a 
and reflection from the region R ~ £, followed by n reflec- 
tions from both the regions R ~ a and R ~ £, followed 
by a final reflection from the region R ~ a and trans- 
mission through the region R ~ £ to the state X. With 
each reflection from the region R ~ £, the amplitude de- 
creases by a factor e~ 2n ' due to the loss of probability 
into high-energy atom-molecule scattering states. 

To calculate the dimer relaxation rate into a specific 
high-energy atom- molecule scattering state X, we need 
the S-matrix element Sx,ad hi ©■ It is extremely sen- 
sitive to the details of the interaction potential at short 
distances through the factor txi- However, the inclu- 
sive dimer relaxation rate summed over all high-energy 
atom-molecule scattering states X is much less sensitive 
to short distances. The matrix elements txi enter only 
in a combination that, by the unitarity of t and the 
explicit expression for t\\ in JSJ, depends only on 77*: 
^2x \txi | 2 = 1 — e~ 4n " ■ Squaring the S-matrix element 
© and summing over X , we obtain 

J2 x \ S X-Ad\ 2 = (1 - e^)| Sl2 | 2 |l - tusu\- 2 . (7) 

This radial law shows that the inclusive dimer relaxation 
rate depends on the interactions only through the param- 
eters a, A*, and 77*. 

To calculate the dimer relaxation rate constant /3, we 
need the behavior of the matrix s as k — > 0, where k is the 
wave number of the atom or dimer in the center-of-mass 
frame. For small k, the most general form for the entries 
of s that is allowed by unitarity and the analyticity of 
Sad, ad as a function of k 2 is 8] 

s n -> -e 2 ^' [1- 2(6 + ib 2 )ak + ...], 

sis -> e^' ' V '4b ak[l - (26 3 + ib x + ib 2 )ak + ...], 

s 22 -> 1 - 2(6 + ih)ak + . . . , (8) 
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FIG. 2: The dimer relaxation rate constant in units of ha/m 
as a function of aA» for r\ t — 0.1, 0.5, 1.0, and oo. 



where bo, b\, bi, and <fi' are real constants. Insert- 
ing these expressions for Sjj into J3J, and expanding to 
first order in k, we find that the (complex) atom-dimcr 
scattering length has the form 



o>ad = (bi — b Q tan[so ln(aA») + <f> + irj*])a . 



(9) 



For 77* = 0, Im a ad vanishes and Eq. JjJJ reduces to Efi- 
mov's expression for the atom-dimer scattering length in 
the absence of deep molecules |8j. The numerical con- 
stants bi — 1.46, bg — 2.15, and <j> = 0.09 were first 
calculated by Bedaque, Hammer, and van Kolck using 
an effective field theory for short-range interactions 

The low-energy dimer relaxation rate can be obtained 
by inserting the expansions JSJ) into the radial law J7J). 
The rate constant j3 defined in is given by the same 
numerical constants bo and <f> that appear in JjJJ: 







20.3 sinh(2?7, 



ft a 



sin 2 [s ln(aA») - 1.48] + sinh 2 77, m 



(10) 



Alternatively, we can use unitarity to obtain (3 directly 
from ©: (3 — — (6irh/m) Ima^D- In Fig. [2 we show the 
dimer relaxation rate constant (3 as a function of aA* for 
a range of 77* . For small 77* , (3 is strongly enhanced at the 
values of a corresponding to an Efimov resonance near 
the atom-dimer threshold. At the same values of a, there 
is resonant enhancement of the threshold cross section 
47r|aAi)| 2 for atom-dimer scattering. As 77* increases, the 
resonance structure is washed out and the coefficient of 
ha/m in i|10|) approaches a constant: (3 ~ 40.5 ha/m . 

For completeness, we also give the universal expres- 
sions for the 3-body recombination rate constants when 
there are deep molecules. Denoting the event rate con- 
stant by a, the loss rate of low-energy atoms from a cold 
gas is dnA/dt = — ian\ . In a BEC, the coefficient of n\ 
is smaller by a factor of 6 [2(j . If a < 0, 3-body recombi- 
nation can only produce deep molecules, while if a > it 
can also produce the shallow dimer. In Ref. [TJ, the rate 



constant a for a > was calculated as a function of a 
and A* if deep molecules are absent. The effects of deep 
molecules on this contribution to a can be deduced using 
the methods described above. In Ref. [l4|. the contribu- 
tion to a from recombination into deep molecules was 
calculated as a function of a and A* for both signs of a 
and infinitesimal 77*. The methods described above can 
be used to generalize those results to arbitrary 77*. The 
two contributions to the rate constant for a > are 
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The only contribution to a for a < is 
4590 sinh(277*) 



ha 4 



[so ln(|a|A*) - 1.38] + sinh 77* m 



(12) 



(13) 



If we set 77* = in (|lll) , we recover the expression derived 
in Ref. 12]. If we expand 1121 |13l) to first order in 77,, we 
recover the expressions in Ref. jlj] after the identification 
7/* = 0.2/ii and correcting for a missing factor of 327T 2 in 
Ref. 01- The expression in 0] for ayeep in the case 
a > has unphysical divergences that are artifacts of the 
expansion to first order in 77*. In the general expression 
(I13J1 . the divergences are replaced by resonance peaks. 

There have been recent measurements of 3-body re- 
combination rates for atoms with large scattering length. 
In Ref. the rate constant a = was measured for 
87 Rb atoms in the |1, +1) hyperfine state in the region of 
the Feshbach resonance near 1007 G. It increases by at 
least two orders of magnitude at the resonance. Since the 
off-resonant scattering length is not large, the expressions 
(lll|) 1131) apply only in the resonance region and there 
were not enough data points in this region to observe the 
scaling behavior a ~ a 4 . In Ref. 0, the rate constant 
a was measured for 133 Cs atoms in the |3, +3) hyperfine 
state in the interval 10 G < B < 150 G, which includes 
several Feshbach resonances. Since the off-resonant scat- 
tering length for 133 Cs is large, the expressions IjllH <|13[) 
should apply except near the points where a goes through 
zero. The data of Ref. provides beautiful confirmation 
of the scaling behavior a ~ a 4 . No evidence of oscilla- 
tory dependence of the coefficient of a 4 on ln(|a|) was 
observed, but neither was it excluded. 

Some recent theoretical predictions of 3-body recom- 
bination rates near a Feshbach resonance have given 
reults inconsistent with the scaling behavior a 



Ref. |21j, the scaling behavior ad e 



K In 

x 2 was obtained 
by assuming a two-step process involving resonant pro- 
duction of the shallow dimer followed by dimer relax- 
ation into deep molecules. Such a two-step process re- 
quires a weak-coupling assumption that breaks down in 
the strong-coupling region where a is large. In Ref. [2l| , 
the authors claim that a scales as a 2 for a > and as |a| 3 
for a < in a particular model for a Feshbach resonance. 
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However their approximations were too crude to extract 
the true asymptotic scaling behavior. 

The rate constants a and /3 in H10(l - l|13[l can be used 
in mean-field descriptions of condensates of atoms and 
dimers near a Feshbach resonance to take into account 
the loss of atoms through 3-body recombination and 
dimer relaxation. While the scattering length a diverges 
as the magnetic field B passes through the Feshbach res- 
onance, the 3-body parameters A* and 77* are smooth 
functions of B that can be approximated by constants if 



the resonance is sufficiently narrow [1J. The study of loss 
processes from condensates could be used to determine 
the parameters A* and 77* for a given Feshbach resonance. 
Once these parameters are known, the binding energies 
and widths of the Efimov states can be predicted and one 
can develop a strategy for discovering these remarkable 
states. 
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